Andreev reflection in bilayer graphene 
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FIG. 1: The graphene bilayer described by Eq. (JT|. 
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We consider the differential conductance of a normal-superconductor junction in clean bilayer 
graphene in the framework of the Dirac-Bogoliubov-de Gennes equation. A remarkable suppression 
of the differential conductance at voltages just below the superconducting gap is found. This can 
be understood in terms of the spinor structures of the electron and hole excitations, in particular, 
the reflected valence-band hole being orthogonal to the incoming electron at normal incidence. 

PACS numbers: 74.45.+C, 73.23.-b, 74.50.+r, 74.78.Na 



I. INTRODUCTION 

Graphene, with its Dirac-like low-lying excitations, 
displays a number of unusual transport properties^! 3 -. 
These result from the pseudospin 1/2 of the low- lying 
modes, their linear dispersion, and the vanishing den- 
sity of states at the Dirac point. The peculiarity that 
graphene retains a finite conductivity down to the Dirac 
p oint^S^^Si ' 10111 allows us to explore unusual regimes. 
Recently normal-superconductor interfaces in graphene 
have been studie d 12 i 13 ' 14 i 15 ' 16 i 17 ' 18 , where the possibility 
of the Fermi energy being smaller than the superconduct- 
ing gap leads to a number of new phenomena based on 
specular— Andreev reflection^. 

In contrast to the monolayer, the low-lying excita- 
tions of bilayer graphene are gapless massive chiral 
fermions with pseudospin L 20 i 21 i 22 i 23 This does not lead 
to such dramatic effects as the absence of backscatter- 
ing in the monolayer ) 24 i 25 i 26 ' 27 which results, e.g., in 
Klein tunneling 2 ^ or perfect electron-hole conversion at 
a normal-superconductor (NS) interface at normal inci- 
dence and all energies below the superconducting gap 12 . 
On the contrary, the pseudospin 1 of the bilayer modes 
is often difficult to connect with observable quantities. 

The purpose of this paper is to examine Andreev re- 
flection at an NS interface in the bilayer and to point out 
observable effects of the pseudospin 1 of the low-energy 
excitations in a way to complement the recently foundi 2 - 
results for the monolayer. As will be pointed out below, 
one key feature is the interplay of the valence-band hole 
mode, which for a given direction has a (two-component) 
spinor structure orthogonal to the electron mode, with 
the other contributing modes. 



II. MODEL 

We consider a clean infinitely wide interface between 
a normal (N, filling the half plane x > 0) and supercon- 
ducting (S, filling x < 0) region of bilayer graphene. The 
bilayer consists of two graphene layers with incquivalent 
sites A and B in one layer and A and B in the other 
layer, with AB stacking (see Fig. [IJ. In the following we 
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with the momentum components k = (k,q), where k is 
normal to the interface and q is parallel to the interface, 
the hopping element t between (nearest-neighbor) A and 
B sites, the chemical potential U, and the Fermi energy of 
the normal region Ep. The factor hv, with v the velocity 
of the monolayer modes, will be set to 1 in the following. 

As will be discussed below, the four-component Hamil- 
tonian {I) is needed to correctly describe the coupling to 
the strongly doped region^ 1 -. This is in contrast to most 
other situations, where the second-order two-component 
Hamiltonian 
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taking into account only the wavefunction amplitudes on 
sites A and B, can be used ^ 20 ' 21 i 22 ' 23 i 28 The advantage 
of the two-component Hamiltonian is the more direct in- 
sight into the pseudospin and the Berry phase 27r of the 
low-energy excitations. This picture may (and will) also 
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be employed for more insight in some particular situa- 
tions in the following calculations, in which the B and A 
amplitudes can be disregarded. 

It is assumed that superconductivity is induced in 
the region S by covering it with a superconducting 
electrode,— or that it is caused by some other intrinsic 
mechanism ) 31 ! 32 so that issues related to the matching 
of the hexagonal lattice in N to a different lattice in S 
(see Ref. |H) do not arise. Further, Eq. (fTJ) assumes that 
there is no layer asymmetry 2 - and also does not consider 
corrections to the dispersion relation which become im- 
portant at very small or large energies^ 

In the absence of a magnetic field, we have time- 
reversal symmetry, and the electron and hole excitations 
of the system are described by the Dirac-Bogoliubov-de 
Gennes equation^ 



H - E F A 
A* Ep — H 
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where u and v are the electron and hole wave functions, 
each having four components denoting the amplitudes on 
the four different atoms in a bilayer unit cell, 
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and e > the excitation energy above the Fermi energy 
E F > 0. 

We assume the normal region to be weakly doped, 
and the superconducting region to be strongly doped, 
described by a potential step of height Uq, 
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(5) 



For large positive Uq, the Fermi wavelength in S is much 
smaller than in N, and the pair potential reaches its bulk 
value at a distance from the interface which is small com- 
pared to the Fermi wavelength in N. This justifies the 
step-function model 



A , x < 
, x > , 



(6) 



where we may choose Ao > real for convenience. The 
junction is assumed to be smooth on the scale of the 
lattice spacing; therefore, the pair potential does not dis- 
tinguish between the sublattices or the layers. In the 
following we will assume the energy scales to satisfy 



U > t > E F , e , A 



(7) 



(the relations between Ep, e, and Ao are not restricted). 
Due to the hierarchy ([7]) it is essential^! to use the first- 
order four-component Hamiltonian ([I]) instead of the 
second-order two-component one @, which is suitable 
for most problems. To obtain a result which is indepen- 
dent of the details of the modeling of the contacts, we 



FIG. 2: The band structure of weakly doped (Ef <C t) normal 
bilayer graphene. At energies below the interlayer hopping t, 
Eq. ((8} gives only one electron (solid line) and one propa- 
gating hole solution. At excitation energies e smaller than 
the Fermi energy E F , the conduction-band hole (dashed line) 
propagates; at excitation energies above the Fermi energy, the 
valence-band hole (dot-dashed line) propagates. 



perform the limit Uq — > oo. The Hamiltonian ([2]), how- 
ever, corresponds to the limit t — > oo, which may not be 
taken first. Therefore the part Uo ^> t of the hierarchy 
(O results in the necessity to use the four-component 
Hamiltonian |T|). 



III. BAND STRUCTURE AND KINEMATICS 

The band structure of bilayer graphene as described 
by the Hamiltonian |1| is given by 



E F = ± t -±\\li 2 -4/,--' 



(8) 



and consists of two vertex-touching hyperbolae accompa- 
nied by vertically displaced high-energy ones. 

In the following we will use the hierarchy of energies 
0, which means that in N only the low-energy bands 
contribute, given the excitation spectrum shown in Fig. [21 
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consisting of a propagating electron mode, the 
conduction-band (CB) hole at e < Ep, and the valence- 
band (VB) hole at e > Ep. The valence band also con- 
tributes an evanescent electron mode, which in general is 
needed to satisfy the continuity conditions at the inter- 
face. For the explicit spinor structures of the modes see 
Appendix [XJ In the following we approximate the spec- 
trum in N as parabolic, e, Ep <C t. The momentum com- 
ponents of the (backscattered) electron are then given 
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FIG. 3: The band structure of strongly doped (Uo S> i) su- 
perconducting bilayer graphene according to Eq. <(3j is given 
by the linear parts of the high- and low-energy dispersion hy- 
perbolae (O, with an excitation gap opened up by the pair 
potential Aq. 
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by 

9 = (e + sin a , (10) 



fc= (c + E F ) cos a. (11) 

At e < (e > Ep), the CB (VB) hole moves opposite 
(parallel) to its momentum, while the respective other 
hole solution is always evanescent. The CB (VB) hole is 
outscattered at the angle 

a = sign (e — Ep) arcsin— j . (12) 

y/\e-E F \ 

In the limits Ep Ao or Ep <C Ao, this leads to retro- 
and specular Andreev reflection, respectively, just as in 
the monolayer situation^ If a > a c , with the critical 
angle a c given by 



a c = arcsin J — - , (13) 
y e + £jp 

the hole cannot be matched to the electron at the in- 
terface using a real momentum. Instead, a' acquires an 
imaginary part and both holes become evanescent. 

In contrast to the normal region, the superconducting 
region is assumed to be highly doped, resulting in the 
excitation spectrum shown in Fig. [3] which results from 
the linear parts of all bands. In the limit Uq — > oo the 
modes in S are normal to the interface. 

IV. NS INTERFACE 

The reflection coefficients of the NS interface can be 
found by elementary mode matching demanding continu- 
ity of the wave function, with the boundary conditions 



FIG. 4: Due to the separation of length scales related to scat- 
tering by the chemical potential and by the pair potential, the 
weakly-doped-N-strongly-doped-S interface may be thought 
of as being decomposed into an NN interface with the poten- 
tial step and a strongly doped NS interface without potential 
step. Scattering processes such as the one shown here are 
summed up to Eq. (|17[) . Solid lines correspond to electron, 
and dashed lines to hole excitations. The scattering matrices 
are given in Appendix iBl 

that the modes in S decay exponentially or propagate 
away from the interface. Equivalently, the effect of the 
interface may be described by a boundary condition^ 3 - on 
the modes in N, which for a strongly doped superconduc- 
tor acts separately on the two layers (since the modes in 
strongly doped bilayer graphene are weakly - symmetri- 
cally or antisymmetrically - coupled monolayer modes), 

v(r) = M u(r) (14) 

with 

/exp{-i/3n.a} \ 

V exp {— ifm ■ a) ) v ' 

and n the unit vector normal to the boundary pointing 
from N to S. Eq. (fT4")) may be used when all excitations 
in S decay (e < Ao) or move (e > Ao) away from the 
boundary. Eq. (|14[) is then valid everywhere in S, so by 
continuity it becomes a boundary condition for the modes 
in N. 

Another way of obtaining the same result arises from 
the observation that the scattering at the junction can 
be separated in two effects happening at well separated 
length scales ; 16 ' 35 The scattering by the chemical poten- 
tial U happens on the scale of the Fermi wavelength of 
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FIG. 5: Reflection probabilities for normal and Andreev 
reflections in the generic situation. Upper graph: Reflec- 
tion probabilities \r\ 2 (electron, solid line) and |»"5 B | (CB 
hole, dashed line) below Ef, plotted against a for t = 10, 
Ef = 0.8 Ao, and e = 0.2 Ao. Lower graph: Reflection proba- 
bilities |r| 2 (electron, solid line) and |f\X B | (VB hole, dashed 
line) above Ef, plotted against a for t = 10, Ef = 0.2 Ao, 
and e = 0.8 Ao. The probabilities of the propagating modes 
add up to 1 (beyond the critical angle a c , the hole is evanes- 
cent). 



the strongly doped region, while the scattering due to 
the pair potential A happens on the scale of the su- 
perconducting coherence length, which in our model is 
much longer than the Fermi wavelength. Therefore, the 
model described above is equivalent to contacting the 
weakly doped normal region by a strongly doped normal 
region, which in turn is in contact to the strongly doped 
superconductor (see Fig. 0|. The NN interface is then 
described by the scattering matrices So(e) for electrons 
and S'o(e) for holes, containing reflection and transmis- 
sion matrices r, r,t,andt, respectively. These matrices 
connect the single propagating mode in the weakly doped 



region (the parabolic part of the lower one of the disper- 
sion hyperbolae) and the evanescent mode existing in the 
weakly doped region with the two propagating modes in 
the strongly doped region (the linear parts of the two 
hyperbolae). It is important to note that at e > Ep the 
relation So(e) = Sq(— e) connects the valence band hole 
scattering matrix with the valence-band electron scatter- 
ing matrix, not with the scattering matrix of the incom- 
ing conduction-band electron. The 2-spinor structure of 
the valence band hole is for a given direction orthogonal 
to the one of the propagating electron. The expressions 
for S and S are given in Appendix [B] for the (more inter- 
esting) case e > Ep, when the reflected hole is from the 
valence band. 

For strong doping, the scattering at the NS interface 
happens between modes which are normal to the inter- 
face. Since Uq ^> t, these modes are weakly coupled 
monolayer excitations. As a result, the subgap scattering 
matrix of the NS interface i o 12 i 16 i 34 



S A (e) = ( 



1 

1 



C = exp 



-i arccos-r- > (16) 
A ' 



in electron-hole space. Summing up multiple Andreev 
reflections between the interfaces in the usual way 3 ^ as 
shown in Fig. 2] results in 



r A (e,a)=Ct[l-C 2 



r r 



(17) 



Note that this is a scalar coefficient since r is a 2 x 2 
matrix and t is a 2 x 1 matrix (see Appendix IB"]). Car- 
rying out the calculations, the reflection coefficients r 
for normal reflection and Ta, corresponding to Andreev 
reflection with the propagating hole solution, satisfy 
\r\ 2 + \i'a\ 2 = 1 at energies below the gap and are plotted 
in Figs. [5] and El 

The differential conductance of the NS junction can 
be calculated from the reflection coefficients using the 
Blonder-Tinkham-Klapwijk formula 36 



§ = 90 (V) 



it/2 



da cosa 



1- \r(eV,a)\ 2 + \r A (eV,a)\' 



(18) 

where go is the ballistic differential conductance of a clean 
sheet of normal bilayer graphene and the reflection coef- 
ficient corresponds to the propagating hole solution in 
N, which, depending on the excitation energy, may be 
from the conduction band (CB) or from the valence band 
(VB), 




CB e < E F , a < a c 
e > Ep , a < a c 
a > a r . 



(19) 



The results are shown in Figs. [7| and [3 The subgap 
conductance shares with the monolayer result the drop 
to zero et eV — A , which just results from the absence of 
propagating holes at any angle of incidence. In addition, 
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FIG. 6: Reflection probability for parameters including a real 
a* as given by Eq. (|2U[1 . Upper graph: Plot of the nor- 
mal reflection probability |r| 2 against the excitation energy at 
the angle of incidence a — 7r/4 for Ef = and t — 3. Before 
falling off to zero (because at eV — Ao right-angle reflection 
is forbidden), jr| 2 becomes unity when a* equals 7r/4. Lower 
graph: Plot of |r| 2 against the angle of incidence a for Ef = 0, 
t — 3, and e = 0.91 Ao in the special interval given by (121[) . In 
this interval, |r| 2 becomes unity at one more angle «*, given 
by Eq. J20J, other than and tt/2. 

a second dip to almost zero occurs in a voltage range 
just below the gap. This extra feature is the main result 
of this paper and is discussed in detail in the following 
section. 



1.5 - 




FIG. 7: Plot of the differential conductance of the NS junc- 
tion, normalized to the ballistic conductance go, for Ef = 
and t = 10Ao (solid line), Ef = 0.5 Ao (dashed line), and 
E F = WA (dotted line). While the dip at eV = E F is 
straightforward to understand, being due to the absence of 
available holes, the dip in the range given by Eq. (I2ip is an 
effect of the pseudospin 1 and destructive interference of the 
two bands, as discussed in the text. 



< e < Ep but with unit probability at e > Ep. The 
physical picture for this is the following: The two-spinor 
(AS) structures of the incoming and the reflected elec- 
tron, as well as of the reflected CB hole, are perfectly 
compatible (the pseudospin 1 does not forbid backscat- 
tering). Andreev reflection at e < Ep is therefore a cor- 
rection of order e/t to this result, arising from the AB 
components. The reflected VB hole, on the other hand, 
has a 2-spinor structure orthogonal to the incoming elec- 
tron. As a result, there is perfect normal backscattering 
at e > Ep. This is an effect similar to the perfect An- 
dreev reflection in the monolayer— (there, both CB and 
VB holes have a spinor structure compatible to the in- 
coming electron, but normal reflection of the electron is 
forbidden). This effect alone fails to translate directly 
into a striking feature in the differential conductance, 
since the faster angular dependence of the spinor struc- 
ture of the bilayer modes washes out those features upon 
angular integration. Under certain conditions, however, 
perfect normal reflection occurs at yet another angle be- 
tween and 7r/2. 

This happens for Ep = at the angle a* with 



V. DISCUSSION 

Unlike the monolayer case, in the bilayer subgap An- 
dreev reflection is always an effect small in e/t. Among 
the characteristic features of the reflection coefficients 
is that, at normal incidence, normal reflection always 
happens with less that unit probability at energies with 



sina* 



tan 2 /3 



1 , (3 = arccos(e/A 



(20) 



for example, at a* — n/4 for t = V3A and e = |v3 A 
(then (3 = 7r/6). The energy interval in which a* is real 
is given by a third-order equation in e; for Ep = and 
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FIG. 8: Plot of the differential conductance of the NS junc- 
tion, normalized to the ballistic conductance go, for Ef = 
and t = 3 An (dashed line) and t = 10 An (solid line). The 
characteristic dip in the range given by Eq. (|2ip moves 
closer to the gap with increasing t. At eV = An, as seen 
from Eq. (IB18|l . g/go = 16/15 (independent of t), while at 
eV < An, the pointwise limit of the conductance is zero as 
t — > oo. 



Aq -C t, it lies close to An. and is given by 



An 6 An 

1 - < — < 1 

2t - A ~ it 



(21) 



where a* = at the lower end of this interval and 
a* = 7r/2 at the upper end. The complete absence of 
Andreev reflection at this angle can be understood in 
the picture of two separate interfaces as the two modes 
available in the strongly doped normal region interfering 
destructively. More insight is obtained by noting that at 
a* , the boundary condition Q14p forms a symmetric dimcr 
on the AB sites if only the CB hole is present (r\ B = 0) . 
This boundary condition is satisfied by r — — 1 and the 
coefficient of the evanescent VB electron mode arising 
from the condition on the A and B sites, resulting in per- 
fect normal reflection at the angle a*. The energy range 
(|2ip thus is the range in which the boundary condition 
(fT"4")l can nontrivially combine the A and B with the AB 
components of the 4-spinor excitations: at smaller ener- 
gies, the orders of magnitude are too different to allow for 
a nontrivial result, leading just to small corrections to the 
e/t — > limit. On the other hand, at energies above (|2"Tj) . 
the mixing angle f3 (which mixes the components of each 
of the two layers but does not mix between the layers) 
becomes too small and the reflection probabilities mirror 
the physics of a pseudospin 1 particle scattered from a 
scalar. 3 - In particular, as can be seen from Eq. (|B18[) . at 
e = A and a — n/i, Andreev reflection becomes perfect 
since scattering off a scalar by it/ 2 is forbidden. The scale 
t thus drops out of the reflection coefficients at e = Ao, so 
that the conductance at eV = Ao depends only on Ep. 
At E F = 0, g/g = 16/15. 



Since in the energy interval (f2"Tj) the probability of nor- 
mal reflection reaches unity at three different values in 
the interval a e [0, n/2], it deviates very little from unity 
over the entire range of angles, leading even after angu- 
lar integration to the characteristic dip of the differential 
conductance shown in Figs. [7] and [5J 

The relevance of this result is that it can be linked to 
the pseudospin of the bilayer excitations, especially the 
forbidden Andreev reflection into the VB hole at normal 
incidence, which, due to the finer angular structure of the 
pseudospin 1 versus pseudospin 1/2, does not drastically 
influence the conductance outside the energy range (|21[) . 

Experimental observations of superconductivity and 
Andreev reflections in monolayer graphene contacted by 
superconducting electrodes have now been reported. 37 
The results given here may therefore be expected to be 
tested in the near future. For the model ((T|) to be appli- 
cable, however, superconducting gaps are needed which 
are larger than the scale below which the Fermi surface 
breaks up into pockets^ The main results presented here 
do not rely on any closer vicinity to the Dirac point other 
than the requirement that the subgap conductance of the 
junction is considered. Therefore, they can be expected 
to be qualitatively valid when there is a connected Fermi 
surface in a part of the subgap energy range, which re- 
quires the superconducting gap Ao, the AB hopping t, 
and the velocities V3 and v describing the AB and direct 
AB hopping, respectively^ to satisfy 3 ^ 



(22) 



Since, compared to the monolayer, the description of the 
bilayer is somewhat more idealized, it would be interest- 
ing to see how the calculated effect survives in experi- 
ment. 



VI. SUMMARY 

In this paper the properties of an NS junction in bi- 
layer graphene have been investigated. The interplay of 
the pseudospin 1 excitations of the weakly doped normal 
region with the strongly-doped superconducting region 
(which may be regarded as two weakly coupled mono- 
layers) produces a nontrivial IV curve. While sharing 
with the monolayer— the features arising from the sit- 
uation Ep < Ao, the most remarkable result is that at 
voltages just below the gap the differential conductance 
is strongly suppressed (see Figs. [7] and [5]). This can, al- 
though indirectly, be traced back to the pseudospin of 
the low-energy excitations in the bilayer. 
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APPENDIX A: BASIS STATES 

In this appendix, we give for reference the electron and 
hole spinors in the normal region. To carry the same cur- 
rent in x direction, these states should be normalized by 
the square root of the x components of their group ve- 
locities which can be obtained from Eq. (jTTJ) . The square 



roots in these expressions follow the convention that the 
branch cut lies along the negative real axis and \J — 1 = i. 
The incoming (left-moving) CB electron is given by 



CB,in 



/ exp{ia} \ 
\ — exp{— ia} J 



(Al) 



There are four outscattered modes in N, the reflected CB 
electron satisfying k 2 + q 2 = t(e + Ep), 



CB, back 



' — exp{— ia} \ 
-^ Jt + E F 

^€ + E F 

i exp{ia} J 



(k-iq) 

- Ve + E F 

(k+iq) 



/ »g-y-g 2 +( e +"E7) \ 



-y /e + Ep 

Ve + E F 

ig+V '-q 2 + {t+E F ) 



(A2) 



the VB electron satisfying k 2 + q 2 = — t(e + Ep), which is always evanescent in our scattering problem, 



VB 



k — iq 
y/ k 2 +q 2 

- We + E F 
-We + Ep 

k+iq 



I -iq+y/-q 2 -(e+E F ) 

- We + E F 
- We + E F 

iq+y/^q~ 2 -(t+E F ) 



(A3) 



the CB hole satisfying fc 2 + q 2 = t(E F — e), 



/ k-iq \ 

~ VE F - e 
VE F -e 

k+iq 



f iq±y/-q 2 + (E F -e) ^ 
VE F -e 

-^JE F - e 



y /E F -e 

iqT\/-q 2 + {E F -e) 
VE F -e 



(A4) 



where the upper (lower) sign is for the propagating (evanescent) CB hole in order to satisfy the boundary condition 
that the mode propagates to the right or decays away from the interface; and finally, there is the VB hole satisfying 
k 2 +q 2 = t(e-E F ), 



k—iq 
y /k 2 +q 2 

-sJ e-E F 
-Ve - Ep 

k+iq 
\ yjk 2 +q 2 ) 



I -i q+ ^- q * + (e-E F ) \ 
Ve-E F 

-Ve - E F 



\ 



'g+\/-9 2 +(e-.EF) 



(A5) 



/ 



APPENDIX B: SCATTERING MATRICES 

In this appendix, we give the explicit expressions for 
the scattering matrices of the weakly doped to strongly 
doped N and the strongly doped N to strongly doped S 
interfaces for energies e with Ep < e < Aq. The propa- 
gating hole is thus the valence band hole, which is the 
time reverse of a valence-band electron. In the limit of a 



large potential step Uq —> oo, the modes in the strongly 
doped region propagate perpendicular to the interface. 
For ease of notation, t and E F are not written explicitly. 



Performing the calculations by elementary mode 
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matching, we find the CB electron-scattering matrix 



So(e) 



m t 

t T r 



m t 12 t 



= I * 



21 ' 22 



12 

22 



(Bl) 



t 21 '22 '22 

[the index 1 (2) denoting the weakly (strongly) doped 



normal region according to Fig. [4] and the index + (— ) 
denoting the upper (lower) band in the strongly doped 
region] with components 



e — yfe cosa — i cos [a — i arcsinh sina] + %\l e(l — sin a) 

m = ■ ■ v (B2) 

— e — v^cosa — zcos [a + i arcsinh sina] — i\J e(l — sin 2 a) 

e 1 / 4 -\/cosa ^— z — i exp{2 arcsinh sina} — expjarcsinh sina}^ 

^21 = ^12 = 7 , \ (B3) 

exp{ arcsinh sina} I e + y^cosa + i cos [a + i arcsinh sina] + i\/e(l + sin 2 a) 



e — -v/e cosa — i cos [a + i arcsinh sina] + i\ e(l + sin a) 

r»~ = V , (B4) 

e + y^e cosa + i cos [a + t arcsinh sina] e(l + sin 2 a) 



r 22 — r 22 + — 



—liyfe sina 



e + ^/ecosa + i cos [a + i arcsinh sina] + i\/e(l + sin a) 



, _h _ 2ie 1 /Vcos5sina 

'12 - *21 - / / v 

e + sin 2 a + i\/e(l + sin 2 a) + cosa ( yfe. + i v 1 + sin 2 a ) 



' 22 



(B5) 



e + -v/ecosa — i cos [a + iarcsinh sinal — i\ e(l + sin a) 

r 2 + 2 + = ; ; V , ■ (B7) 

e + Y^cosa + i cos [a + zarcsinh sina] + i\ /e(l + sin 2 a) 



I 

Direct inspection verifies that S'o(e) is symmetric and uni- band. We find 
tary. 

As mentioned in the main text, the VB hole scattering _ / fii i 12 *ia 

matrix is related to the valence-band electron-scattering -So(e) = t 21 f 22 r 22 + (B8) 

matrix by So(e) = «Sq(— e) and thus needs to be calcu- \ t 21 r 22 

lated separately since the spinor structure of the valence 

band differs from the spinor structure of the conduction with components 
I 



e — y/e cosa — i cos [a — i arcsinh sina] + i\/ e(l + sin a) 

hi = V . (B9) 

— e — -y/ecosa — zcos [a + i arcsinh sina] — i\ /e(l + sin 2 a) 



— 2ie 1 / 4 v / cosa sina _ , 

*ak = *ia = / ■ V ~ ( : r (BIO) 

e + sin 2 a + i\Je(l + sin 2 a) + cosa ( y/e + i\J 1 + sin 2 a I 



e + ^/ecosa — i cos [a + i arcsinh sina] — iy e(l + sin a) 



e + \/~€ cosa + i cos [a + i arcsinh sina] + i \ / e(l + sin a) 



(Bll) 



-2iyf£ sina 



'22 — ' 22 



e + -y/e cosa + i cos [a + i arcsinh sina] + iy e(l + sin a 



(B12) 



f+ _ f+ 



e^-y/cosal i + i exp{ 2 arcsinh sina} + 2y / eexp{arcsinhsina} 



exp{ arcsinh sina} e + y/e cosa + i cos [a + z arcsinh sina] + iy e(~l + sin a 



(B13) 



e — ^/ecosa — i cos [a + i arcsinh sina] + iy e(l + sin a 



e + ^/ecosa + i cos [a + i arcsinh sina] + iy 1 e(l + sin 2 a) 

r 



(B14) 



As Sq, also So(e) is symmetric and unitary. Inserting the 
above expressions into Eq. (fTT|) gives the results plotted 
in the figures of the main text. 

While the corresponding analytical expressions for r 
and ta in the general case are too lengthy to be displayed 
here, the following limiting cases can be given: 

(1) For Ep = 0, leading order in e/t <C 1 (the scale Ao 
drops out): 



r = 1 — 2 sin 4 a + 2i sin 2 avl 



• 4 

sin a 



rT = 0, 
rT = 0. 



= 1. 
(B15) 
(B16) 
(B17) 



(2) For E F = 0, e = A [this means (3 = and M = I, 
now the scale t drops out since M does not mix the large 
(AB) and small (AB) components of the 4-spinors]: 



-icos2a 



2 / • 4 

-i sin a + v 1 — sin a 



-1 + sin 2 a — iyl 

— z sin 2 a + yl 
2i cosa sina 



• 4 

sin a 



sin 2 a 



• 4 

sin a 



|r| 2 = cos 2 2a , (B18) 

- , (B19) 
(B20) 
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